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Abstract. CO J = 5 — 4 and J = 2 — 1 emission lines 
were detected towards the extremely red galaxy (ERG) 
HRIO (J164502-h4626.4) &t z = 1.44. The CO intensities 
imply a molecular gas mass M(H2) of 1.6 x 10^^ h^Q^ M©, 
and, combined with the intensity of the dust continuum, a 
gas-to-dust mass ratio around 200-400 (assuming galactic 
values for the conversion factors). The peak of the CO 
lines are at the same redshift as the [On] 3727 line, but 
blue-shifted by 430 km s~^ from the Ha line. These CO 
detections confirm the previous results that HRIO is a 
highly obscured object with a large thermal far-infrared 
luminosity and a high star-formation rate. The overall 
properties of HRIO (CO detection, LpiR to L^q ratio, and 
FIR to radio flux ratio) clearly favour the hypothesis that 
its extreme characteristics are related to star-formation 
processes rather than to a hidden AGN. 
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1. Introduction 

The recent detections of CO emission at cosmological dis- 
tances provide hints about the physical structure of newly 
formed objects (Combes, Maoli & Omont 199£ and refer- 
ence therein). A measure of the total gas and dust mass is 
indeed a very useful indicator of the object evolutionary 
status, because it provides an estimation of the fraction 
of the galaxy which has yet to be turned into stars at 
the epoch of observation. At high redshifts such measure- 
ments, therefore, provide hints about the occurrence of 
active star-formation processes and help in investigating 
models of galaxy formation (Silk & Spaans 1997). 

Only a handful of distant objects was detected so far 
in CO, and most of them appear to be magnified by grav- 
itational lenses (table 1). For these objects the mass of 
molecular gas inferred from the CO intensities - corrected 
for gravitational amplification and using the Galactic CO 



to H2 conversion factor - turns out to be 30 80 % of 
the total dynamical mass, with typical values of M(H2) 
~ 10^° 10^^ Mq. These masses could be somewhat 
smaller if the CO to H2 conversion factor were higher than 
in the Galaxy. Large quantities of CO are expected when 
large FIR luminosities and dust content are detected. Most 
of the objects detected until now in CO were indeed se- 
lected because they were dust-rich systems. 

HRIO belongs to the class of objects with very red 
colours {R~K > 6). Their faintness at optical/NIR wave- 
lengths makes the redshift determination and the investi- 
gation of their nature difficult even with 4m class tele- 
scopes. HRIO is so far the only one with a measured spec- 
troscopic redshift (Graham fc Dey 1996). One of the main 



issues regarding ERGs is whether they are young and star- 
bursting galaxies hidden in the optical by a large amount 
of dust or whether they are old passively evolving galaxies 
at z > 1. Recent results indicate that both classes con- 



tribute to the population of ERGs ( Cimatti et al 1999 ). 
HRIO was first detected in the submm/mm continuum 
with the IRAM 30m equipped with the MPIfR bolome- 
ter and with the JCMT equipp ed with the SCUBA 
double arrays ( Cimatti et al 1998 ). Subsequent observa- 
tions confirmed the submm detection of this galaxy 
(Dey et al. 1999). The inferred properties of this object 
show its extraordinary nature: its dust mass is 4 8 x 10® 
Mq and its total FIR luminosity in the range 10 -2000 
^m rest-frame is 2 2.5 x 10^^ Lq (Hq = 50 Mpc/km/s). 
This places HRIO in the class of the ultra-luminous in- 
frared galaxies and suggests the presence of a star-forming 
object with a SFR of - 200^500 Mq yr'^ or even higher. 

In this paper we present C0(2-l) and CO(5-4) obser- 
vations of HRIO made with the IRAM Plateau de Bure 
interferometer. The observations are described in Sect. 2, 
the resulting detection of both lines are presented in Sect. 
3, while implications of these measurements are reported 
in Sect. 4. Throughout the paper, we adopt Hq = 50 
Mpc/km/s and qo = 0.5. 
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Table 1. List of high-redshift objects detected to date in CO. 



Name 


type 


redshift 


(10" Mq) 


references 


















F10214+4724 ' 


ULIRG 


2.28 


4.2* 




Solomon et al. 1992 






Cloverleaf ' 


QSO 


2.558 


0.9* 




Barvainis et al. 1994| 






BR1202-0725 ' 


QSO 


4.69 


2.4 




Ohta et al. 1996, 


Omont et al. 1996 




MG0414+0534 ' 


QSO 


2.639 


3.3 




Barvainis et al. 199^ 






SMM02399-0136 ' 


submm-HLG 


2.808 


1.8* 




Prayer et al. 1995 








SMM14011+0252 ' 


submm-HLG 


2.565 


1.1* 




Prayer et al. 199^ 








APM08279+5255 ' 


BAL quasar 


3.911 


0.03* 




Downes ct al. 1999 






BRI1335-0417 


QSO 


4.407 


2.4 




Guilloteau et al. 1997| 






53W002 


weak RG 


2.39 


1.3 




Scoville et al. 1997a|, |Alloin et al. 1999 




4C 60.07 


RG 


3.79 


1.8 




Papadopoulos ot al. 1999 






6C 1909+722 


RG 


3.53 


1.0 




Papadopoulos et al. 1999 







I: evidence for a lensed object 
*: corrected for magnification 



2. Observations 

The observations were done partly during the winter 1998- 
99 and partly during the summer of 1999 with the IRAM 
Plateau de Bure Interferometer. A dual SIS 3mm/1.3mm 
receiver was used to observe simultaneously the CO (2-1) 
(redshifted in the 3mm band) and the CO (5-4) (redshifted 
in the 1.3mm band) lines. In both cases, the 3mm receiver 
was connected to two units of the Correlator, providing a 
velocity coverage of about 850 km s^^, while 4 units were 
connected to the 1.3mm receiver providing a velocity cov- 
erage of about 650 km s~^. 

The observations made during the winter 1998-99 assumed 
a redshift of 2 = 1.443 deduced from the Ha line. This 
yields frequencies of v = 94.405 GHz for the C0(2-l) hue 
and V = 235.982 GHz for the CO(5-4) line. Both lines 
were detected, but blue-shifted by about 400 km s~^ with 
respect to the Ha line and thus truncated. The observa- 
tions were repeated during the summer 1999 with a new 
frequency setup: ly = 94.531 GHz and v = 236.297 GHz 
for the C0(2-l) and CO(5-4) hues respectively. 
The observations were made in the CD configuration of the 
interferometer. The phase drifts were calibrated by observ- 
ing the nearby quasars 1633-1-382 and 1732-1-389 every 20 
minutes throughout the observations. The amplitude was 
calibrated with the quasars 3C273, 3C345 and 2145-1-067, 
and the compact Hii region MWC349 at the beginning 
and/or the end of each transit. The passband of the sys- 
tem was calibrated using 3C273 or 2145-1-067. At 1.3mm, 
the total useful integration time on source was about 9 
hours during the winter 1998/99, and 12 hours during the 
summer 1999. More integration time was available at 3mm 
because some additional data were obtained, marginal at 
1.3mm but of good quality at 3mm. 

The 3mm receiver was tuned in single side band, while the 
1.3mm was tuned in double side band, with the CO(5-4) 



in the lower side band. The upper side band was then 
available to measure the continuum. 

Three visibility tables were produced from the cali- 
brated data: two spectral tables containing the C0(2-l) 
and CO (5-4) data respectively and one continuum ta- 
ble with the 1.3mm continuum. The spectral data were 
smoothed in frequency to improve their signal-to-noise 
ratio: the velocity resolution of the C0(2-l) data was 
smoothed to 50 km s~^ and that of the CO(5-4) to 75 km 
s~^. The UV tables were imaged, CLEANed and restored 
with eUiptical gaussian beams of 5.6" x 4.6" {PA = 56°) 
and 2.9" x 2.0" {PA = 65°) for the 3mm and 1.3mm data 
respectively. 

3. Results 

Signal is detected in the central channels of both the 
C0(2-l) and CO(5-4) spectral data sets at 10 and 5 sigmas 
respectively. The averages of the channels with detected 
emission were used to produce the integrated maps shown 
on Figs, la and Ic. The C0(2-l) integrated flux is 1.4 Jy 
km s"^. There is no indication for any extension of the 
emission at this resolution. The peak of the source is off- 
set by about 1" from the position of the optical source 



as measured by HST (Dey et al. 1999). This offset is well 



within the astrometric accuracy of HST but if real it may 
be indicative of a spatial extinction. 

The corresponding CO (2-1) profile (Fig lb) is roughly 
gaussian, with a FWHM of almost exactly 400 km s^^, 
and a central frequency that corresponds to a redshift of 
z = 1.439 ± 0.001. The appearance of a flat-top or even 
double peak profile around ~ 150 km s~^ (see Figlb) 
cannot be checked with the present data and only ob- 
servations at higher signal-to-noise ratio (those exploiting 
the full spectral resolution) can settle its reality. The line- 
width observed is fairly large, although not atypical for 
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Fig. 1. (a) C0(2-l) integrated intensity map; the contours are at 0.25 Jy beam~^ km s~^. (b) C0(2-l) spectrum 
at the center of the source, (c) CO(5-4) integrated intensity map; the contours are at 0.75 Jy beam~^ km s~^. (d) 
CO (5-4) spectrum at the center of the source. The contribution from the continuum is not subtracted from the Hne 
flux. 



this kind of sources (e.g. SMMJ02399, Frayer et al. 1998) 
and could be due either to an edge-on system or to various 
separate components. The redshift deduced from the CO 
line is apparently shifted from that deduced from the Ha 
line (z = 1.443), but corresponds to the redshift deduced 
from the [Oil] line at 3727 A (z = 1.439). 
This shift is still within the uncertainties in the opti- 
cal redshift but if real it would be different from what 
i s typically found for l ow redshift luminous galaxies 
( ganders fc Mirabel 1996| ), where systematic blue- ward 
offsets of optical lines from the CO redshift are at- 
tributcd to outflows with d ust obscuration (see e.g., 
Gonzalcz-Dclgado et al. 1998 ). The question remains. 



however, completely open since a recent analysis by Mcin- 



tosh et al. (1999) of a sample of quasars shows how high- 
redshift objects present H/3 lines with a systematic mean 
red- ward shift of ~ 500 km/s with respect to the systemic 
redshift of the objects (that defined by the narrow line re- 
gion). Even though the comparison with quasars may not 
be fair since the line emitting regions could be different, 
HRIO seems to show similar properties with the CO red- 
shift corresponding to that of the narrow forbidden lines 
and coinciding with the centre of mass of the system, while 
the Ha line is shifted with respect to that. 

The integrated map corresponding to the CO (5-4) map 
(Fig. Ic) shows a source at the same position as the C0(2- 
1) integrated map. The integrated intensity of that source 
(corrected for the contribution of the continuum) is 1.35 
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Jy km s~^. The corresponding CO(5-4) profile (Fig. Id) 
is roughly gaussian, with a FWHM of 380 km s~^ similar 
to that of the C0(2-l) profile. The central frequency cor- 
responds to a redshift oi z — 1.440 ± 0.001 similar to that 
deduced from the C0(2-l) fine. 

Although the CO J = 5 level is J( J + 1) • 2.77K 83K 
above the ground state the integrated flux (in Jy km 
s-i) of CO(5-4) is equal to that of the C0(2-l) line, 
(5 — 4)/(2 — 1) ^ 1. CO luminosities, in solar units, are 
1.5 and 3.7 x 10^ L© for the C0(2-l) and CO(5-4) line re- 
spectively, while the total line luminosities L^q are 4 x 10^" 
and 6 X 10^ K km pc^ h^^, for the C0(2-l) and C0(5- 
4) line respectively. When expressed in these latter units, 
the ratio between the C0(2-l) and CO(5-4) luminosities 
of the same source is proportional to the line intrinsic 
brightness (Rayleigh- Jeans) temperature ratio integrated 
over the area of the source: 

rb[CO(5-4)] _ £^o(5-4) a(2-l) 
rb[C0(2-l)] i^o(2-l) ^s{5-4) ^ ^ 

where fig is the source solid angle. If the spatial extent of 
the CO(5-4) emission region is similar to that of the C0(2- 
1) - a plausible hypothesis since both transitions have 
same line-width and profile - = 0.57 and corresponds 
to a value of the excitation temperature of Tox 18 K 
(see e.g. Maloney, &: Black 19881 ). If gas and dust are in 



thermodynamic equilibrium the kinetic temperature Tkin 
would equal Tdust, but it is usually found that Tkin < Tdust- 
I n HRIO the dust temperature was es timated to be ^ 40K 
(ICimatti ct al IQQSj |Dey et al. 199g| ) and this value can 
be taken as the upper limit to the gas temperature. If 
we assume that Tkin — 20 K the gas density implied by 
this ratio is less than 10"^ cm^"^. As an example, at Tj-ad ~ 
2.77(1 + z) = 6.76 K and T^in = 20 K the estimated 
excitation temperature is Tex = 11 K for a CO density 
of 300cm~'^. The present data cannot distinguish between 
a picture where the dominant component of the ISM in 
this system is a diffuse {n{H2) ^ 10^ — lO'^cm"'^) gas or 
whether the medium is clumpy. The spatial shift between 
the CO and Ha lines would be more compatible with this 
latter picture. 

If the ratio between L'^q and the mass of molecular 



gas is similar in HRIO and Arp 220 ( |Scoville et al. 1997bD , 
the molecular gas mass in HRIO, using the C0(2-l) 
line, is M(H2) = 1.6 x 10" Yi^^Mq, larger than what 



is usually found in local ULIRGs (3olomon ct al. 1997 
Braine & Dumkc 1998) but similar to that of other de- 



tected high-z sources (Frayer et al. 1998). 

The 1.3mm continuum map shows only a marginal 
(2ct) detection at the position of the source, with an 
integrated fiux of 2.2 ± 0.9 mJy beam^^. The flux de- 
tected with the IRAM 30m telescope was 4.9 ± 0.8 mJy 



(Cimatti et al 1998) at 240 GHz. Scahng the flux as v'^ 



the PdBI detection would correspond to an expected flux 
of 2.4 ± 1.0 mJy at the observed frequency of the 30m. 
We exclude that this discrepancy is due to an extended 



component; it is more likely due to the difference in the 
calibration of the two instruments since the two values are 
consistent within the error bars. 

4. Implications 

The analysis of the spectral energy distribution of HRIO 
shows the presence of thermal emission at rest-frame 
A > 60^m with a range of dust temperatures be- 
tween 30 and 45 K. The implied total dust mass is 
8 — 4 X IQ^h^n M0 (for a dust emissivity index (3 of 
Therefore the resulting gas-to- 



2, Cimatti et al 199? 



dust mass ratio for HRIO ranges between 200 and 400, 
as local spirals ( Andreani, Casoli fc Gerin 1995| ), ULIRGs 



(Solomon et al. 1997) and also sub- mm selected luminous 
sources show (Frayer et al. 1999). 

The total rest-frame far-IR luminosity in the range 10 — 
2000/iTO is 2 - 2.5 x lO^'^h^^ Lq (pimatti et al 1999D as 



estimated taking into account the ISO upper limits at 90 
and 170 ^.m (Ivison et al 1997). When these latter are not 
considered and the 450 //to detection is i ncluded the lum i- 
nosity turns out to be a factor of 3 larger ( Dey ct al. 1999 ). 
The ratio lies therefore in the range ~ 60 — 175iQ 

(K km s^^ pc^)~^, which agre es with the relation foun d 



for nearby luminous galaxies ( Sanders fc Mirabel 1996 ) 



whose emission is mainly powered by star-formation. Ob- 
jects whose FIR emission is dominated by an AGN - as the 
hyperluminous Infrared Galaxies - show much larger ^T™ 

and do not even show up in CO (e.g., Evans et al. 1998). 



This indicates that the overall FIR emission by HRIO is 
dominated by star formation. Assuming that most of the 
FIR luminosity is due to recent OB star formation ac- 
tivity, the star formation rate turns out to be SFR — 
^ 10-i°Lfir ~ 200 - 500/i^(f Mg/yr. 
Star formation efficiency is usually measured by the ra- 
tios and (see e.g. [Young 1999| ). While the for- 
mer shows indeed quite a high value (16 — 44) similar to 



that of merging local systems (Young 1999), the latter is 
of only 0.007 and very likely indicates a large extinction 
affecting the Ha emission. 

With the values above for molecular mass and SFR this 
active phase of gas depletion lifetime should have lasted 
at least: 



<gas = 4 X 



10^ 



5(:^)-i = (0.2-0.6) 



10« 



(2) 



CO 



The large value of gas conversion into stars (with re- 
spect to local galaxies) could be consistent with two possi- 
ble scenarios: either a genuinely young galaxy in the pro- 
cess of active star-formation (and the detected amount of 
gas seems enough to feed it), or the presence of a large 
amount of gas could be the result of a merging process of 
two discs (in this latter case the resulting galaxy will have 
a mass of a present-day massive elliptical) . 

Most of the properties of HRIO suggest that the 'lo- 
cus', which best characterizes it, is that of local ULIRGs 
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(Hughes, Dunlop, Rawlings 1997). HRIO follows also the 
expected tight correlation between the infrared flux and 
the radio continuum: in fact the logarithmic ratio of FIR 
(60 ^m) and radio (1.5GHz) continuum flux density (in 
HRIO rest-frame) q — log j™ ~ 3 again falls within 
the value of nearby starbursts ( ganders fc Mirabel 1996| ) . 
Furthermore, the ratio between the line (2-1) and (5-4) 
and the FIR luminosities, ^2212^ = (2 - 6) x 10"^ and 

' -Lfir ^ ' 

= (0.5 - 1.5) X 10-^ agree with a model of CO 
emission in high redshift galaxies, based on an extrapola- 
tion of the properties of local ULIRGs ( Blain et al. 1999D . 

With the detected line width and the upper limit on 
the CO source size, given by the effective beam width 6 < 
5", the upper limit to the total dynamical mass contained 
within the CO emitting region is: 

A^dyn = i^f < 7.7 X lO'\sim)-^h^,'MQ (3) 
z KjT sm I 

where Ay is the observed deconvolved line width, i is 
the inclination and R is the linear diameter of the source 
(i? < 43/i^Q^A:pc). The resulting dynamical mass is a factor 
of 5 larger than the estimation of the molecular mass. 
The two values would coincide if the CO emission were 
concentrated within the inner 10 kpc. 
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